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ABSTRACT 



An experimental study has been made of the feasibility of 
following the precipitation process in a polycrystalline aluminum- 
copper alloy. The intermediate structures formed during ageing 
have previously been determined by single crystal methods to pass 
through the following sequence* 

G.PAO — G.p.ui — e' — e 

The sequence of structures was followed by analysis of the diffrac- 
tometer line profile of the ©<.('2.00) peaks. The experimental 
profiles were corrected for instrumental broadening by the Stokes 
method. One-dimensional Patterson projections of the electron 
density distribution in the alloys were synthesized from the 
corrected profiles. 

The complete sequence of intermediate precipitates could not 
be realized in this investigation. However it was possible by this 
method to distingish between the G*P. and Q structures by 

observing the height and shift of the Patterson function origin peaks* 
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1 . 



Introdtic tion. 



The phenomenon of precipitation hardenixig holds great commer- 
cial interest. It is a method of enhancing the mechanical properties 
of a metal without resorting to other processes such as work hardening, 
dispersion hardening or a martensitic type transformation. It elim- 
inates the concern regarding the hardenability of a metal. Metal 
parts may also be formed or fabricated in a soft condition and then 
hardened by a relatively simple ageing treatment while retaining 
excellent dimensional stability. One example of the importance of 
precipitation hardening is its application to aluminum alloys which 
has resulted in a great number of usable high strength materials. 

The mechanism of precipitation hardening is not at all simple. 

The most basic concept involves the formation of a stable phase 
dispersed throughout a matrix. This does not include dispersion 
hardening, which is due to the presence of a second phase dispersed 
in the matrix. However, it is not a precipitation process. Even 
in the basic concept of precipatation there are constraints imposed 
on the nucleation and growth of such a phase because it is being 
formed in a solid solution. In general, the process is much more 
complicated than this. It varies with each alloy system but usually 
involves the precipitation of a metastable phase, which later gives 
way to a stable one. There may be as many as three intermediate 
stages before the stable precipitate is formed. The intermediates 
may or may not be true phases, depending on the alloy system. 

The intermediate stages hold great interest. In many alloy 
systems these intermediates have the greatest influence on mechanical 
properties, and the formation of a stable precipitate results in a 
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decrease in hardness and strength. 

The Al-Cu alloy represents an ideal system in which to study 
the age-hardening process. It is representative of many systems for 
it displays both substitutional and displacement disorder during 
ageing. The copper atoms have a greater scattering factor and a 
smaller radius than the aluminum atoms = 29/l3 and 2.56 A 

compared to 2.86^ )• The difference in atomic size results in the 
displacement disorder while the greater scattering factor of copper 
allows its effects to be detected by x-ray analysis. 

The age-hardening system in Al-Cu has been studied more than 
any other. As a result its sequence of structures is knoivn in consid- 
erable detail.^ 

The mechanism of the decomposition of the quenced supersatur- 
ated solid solution shows four different stages whose relative roles 
depend strongly on the annealing temperature. In general the sequence 
may be represented as follows: 

G.R 0] e' e 

Two stages, 9 and0, correspond to true precipitates of composition 
CuAl^ , where 0 is the equilibrium phase. The G.P. refers to a 
’’Guinier-Preston zone”, which may be defined as a region rich in 
solute atoms which is formed on certain planes in the solvent matrix 
and retains what is essentially the structure of the matrix. G.P.^\^ 
is a copper rich region of plate-like shape formed on planes 

of the aluminum matrix whose maximum size is a few atomic planes 
thick and up to about 100 A in diameter. G.P. displays features 

Ij. M. Silcock, T. J. Heal and H. K. Hardy, J. Inst. Metals 82 , 
240-242 (1953). 
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of a zone and a true phase. As these platelets reach a thickness 

o 

of about 30 A they become unstable and dissolve. They may be considered 
as having a lattice similar to that of the aluminum matrix. To a first 
approximation it is merely a superlattice of altiminum, tetragonal 
with c/a = 2 where a is the parameter of aluminum. 

The above structures do not necessarily evolve from one another 

in sequence, but two (or more) may occur simultaneously in the same 

specimen. The conditions have been determined in which the three 

metastable phases and the equilibrium phase appear as a function of 

o 

time and temperature for different concentrations of Cu in Al. 

A variety of x-ray techniques have been used in the numerous 
precipitation investigations. Single crystals have generally been 
employed using both white and monochromatic radiation. Some work 
has been done using normal powder methods. The small angle scatter 
method has been applied with great success by Guinier.^ The tech- 
niques using single crystals and monochromatic radiation have 
advantages of ease and certainty of interpretation in terms of a 
reciprocal lattice plot, but require sophisticated equipment and 
precise experimental technique. Some of the results obtained by 
the above methods have recently been confirmed by electron micro- 
scopic photographs."^ The use of an x-ray diffractometer in the 
analysis of polycrystalline specimens is attractive because of the 
relatively simple procedures involved and the possible extension of 

^ Ibid . , pp. 243-45. 

Q 

A. Guinier, Ann. Phys. 1^, 161 (l939). 

^G. Thomas, Transmission Electron Microscopy of Metals , pp. 236-248. 
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the technique to precipitation studies in commercial alloys. Hence, 
it is the purpose of this investigation to follow the ageing process 
in a polycrystalline Al-4^Cu alloy by means of a diffraction line 
shape analysis and Patterson projection of the electron density 
distribution within the alloy. 

The data necessary to form a Patterson function is available 
from a diffraction line profile, as is shoim below. However, a true 
line profile must first be obtained by correcting the observed profile 
for instrumental broadening, due mainly to the slight spread of wave- 
lengths in the incident radiation. A method has been derived by 
Stokes using the convolution theorem of mathematics to correct the 
experimental line shape without any prior assumptions concerning the 
corrected line shape. The method consists essentially of obtaining 
the Fourier coefficients for the broadened and reference curves, 
operating on these coefficients (termed unfolding), and then synthe- 
sizing the corrected curve from the unfolded coefficients. This 
technique eliminated broadening due to Kol doubling as well as other 
forms of broadening. 

It has been shoim that the true line profile is proportional to 
the Fourier transform of the square of the absolute magnitude of the 
structure factor, The structure factor general a 

complex quantity, represents something real about a crystalline 
material for it is the Fourier transform of the electron density in 



^A. R. Stokes, Proc. Phys. Soc. London 382 (l948). 
^A. J. C, Wilson, X-Ray Optics , p. 107. 
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the material*^ Experimentally only determined which yields 

no information concerning its phase* If the phase relations could be 
determined experimentally there would be no problem in finding the 
arrangement of atoms in any crystal, no matter how complicated. Of 
the methods proposed to circumvent the phase problem, the Patterson 
function has emerged as the most powerful and adaptable. This function 
is a Fourier series whose coefficients are If one considers 

the Fourier coefficients as derived by forming products of 

each Fourier coefficient Fv\VJL with its complex conjugate, the electron 
density function can be derived directly from the Patterson function; 
The Patterson function is essentially a display of all the vector 
lengths between atoms of a unit cell. All the vectors of one length 
will then appear as a peak in this display. Since the reverse vectors 
are also displayed, it is necessary to sort out the duplicated peaks 
in order to arrive at the electron density function. 



2. Alloy preparation. 

The Al-4^Cu (by weight) alloy was prepared using 99.99^ ingot 
A1 and 99.94^ Cu shot. The alximinum was first melted in an induction 
heated graphite crucible, and then the copper shot was added. The 
molten metal was held between 700® and 720® C for 45 minutes to ensure 
complete solution of the copper. The alloy was cast into a flat 
rectangular zirconia boat, prehee.ted to 700® C, and allowed to solidify 
at room temperature. The ingot thus formed weighed 125 grams, 

7 

A. Guinier, X~Ray Diffraction , p. 92. 

g 

M. J. Buerger, Vector Space , pp. 5-10. 
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measuring approximately 4” x 2 ^'^ x 

Microscopic examination of the ingot revealed a large grained 
structure containing intergranular theta phase and no evidence of 
undissolved copper. 

The ingot was homogenized at 542^ C for 24 hours and furnace 
cooled. Grain refinement was accomplished by cold rolling and 
recrystallizing at 600^ C for one hour. The amount of reduction 
per pass through the rolls was limited to about 20^ due to work 
hardening, so the refinement treatment had to be repeated three 
times, resulting in an ingot .17” thick, and containing equiaxed 
small grains with finely divided theta phase dispersed throughout. 

Eight hardness specimens were cut from the ingot, each being 
about 1.75” X .5” x .17” in size. A hole was drilled in one corner 
of each specimen and top and bottom surfaces were smoothed on no. 0 
and no. 3/0 metal lographic paper. 

3. Solution treatment and age hardening. 

The eight specimens were strung on a wire with a separation 

to ensure rapid heat transfer both on heating and quenching. They 

o 

were solution treated in a salt bath for six days at 536 C, which 

O fff o 

is about 36 above the solid solubility of 4yo copper in aluminxim. 
Following the solution treatment the specimens were quenched in 
room temperature water, and the as quenched hardness was taken 
immediately on one specimen. The surface irregularities caused by 
the salt bath on the remaining specimens were removed by lightly 

%f. Hansen, Constitution of Binary Alloys , p. 84. 
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polishing on no. 0 and no. 3/0 paper. They were then suspended in 
an oil bath at 165^ C for ageing. The as quenched hardness specimen 
was aged at room temperature. 

Microscopic examination of the as quenched hardness specimen 
following solution treatment revealed a small amount of undissolved 
theta phase in the structure. A lattice parameter measurement on 
the oC(422) diffraction peak of this specimen indicated that about 
2^% of the copper had gone into solution. The small amount of undis- 
solved theta would have no effect other than to reduce the peak 
hardness obtained during ageing. 

4. Diffraction procedure. 

As each specimen reached its respective ageing time, it was 
quenched in room temperature water and hardness readings were taken 
immediately. No less than six readings were taken on both sides to 
ensure a good average. The surfaces were then given a rough electro- 
polish to remove the oxide layer and to ensure that the surface was 
as nearly a true representation of the bulk of the specimen as possible. 
Chemical macro-etching was not used due to its preferential removal of 
copper from the alloy. 

The 06 ( 200 ) peak of each specimen was recorded using a Norelco 
x-ray diffractometer with Ni filtered CukoL radiation. The diffracted 
radiation was counted by a xenon proportional tube in conjunction with 
a pulse height analyzer set to pass the central 90^ of the pulse 
distribution. 

The reference peak for the Stokes* unfolding was obtained from 
a recrystallized specimen of the pure ingot aluminum. It was possible 
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to get a well resolved koc doublet for the (200 ) peak from this specimen* 
5* Treatment of data. 

The intensity of each diffraction peak was recorded as a function 
of "2.0, Due to the sharpness of the reference peak it was necessary to 
use an interval of .0125^20 between points to describe the peak by a 
smooth curve. It was possible to use a point interval of .025*^ 2.9 to 
describe the broadened peaks. A total interval of 1.8^ was used on 
all broadened peaks except the as quenched and aged one day specimens, 
where it was possible to use slightly shorter intervals due to their 
narroimess. The rnaximum value of the diffraction peak was centered 
in the total interval in all cases. 

The unfolding and Patterson function synthesis calculations were 
performed by the TT, S. Naval Postgraduate School’s CDC 1604 electronic 
computer. The diffraction intensity data was punched onto computer 
cards for input to a Fortran program developed by Prof. J. R. Clark, 

Dept, of Material Science and Chemistry, to execute the above calcu- 
lations. 

6. Experimental results. 

The hardness/ageing time curve for the specimens aged at 165^ C 
is show in Fig. 1, page 9. The specimen aged at room temperature for 
65 days attained a Rockwell ”F” hardness of 77.4. 

The lattice parameter of the as quenched specimen was 4. 044 A > 
which corresponds to 2i weight per cent copper in soltition. The 
specimen aged at 165 C for one and two days retained this same lattice 
parameter. The lattice parameter of the specimens aged for five days 
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and longer was 4, 049 A , which is the same as that of pure altoninum. 
The as quenched specimen retained the same lattice parameter on ageing 
at room temperature for 65 days. 

The Patterson functions for the artificially and naturally aged 
specimens are shown in Appendix I on both a small and large scale. 

The ordinate is P(w), where is the interatomic vector length. The 
vector between any two atoms in the crystal always has its tail at 
the origin in the Patterson function. The abcissa then represents 
the magnitude of w as it ranges from zero to ”a”, the lattice param- 
eter of alimiinum, in increments of .Ola. The Patterson functions 
displayed in Appendix I are the one-dimensional projections of the 
electron density in the direction. 

7. Discussion of results. 

The lattice parameter and hardness measurements were used to 
confirm the structures present at each ageing time. The lattice 
parameters were determined with only sufficient precision to detect 
any changes as a function of ageing time. The structures contained 
in each specimen were determined by comparing the experimental 
hardness curve with the hardness cxu'ves of Silcock, et al.^^ It 
was not possible to obtain the complete sequence of intermediate 
precipitates on ageing at 165^ C since only copper went into 

solution. Initial ageing at 165^C produces G.P. and © simul- 

taneously, while at peaJi hardness (five days) 9 predominates. The 
specimen fully aged at room temperature contains only G.P. \X\ since 

^^Silcock, et al. , op . cit . , p. 243. 
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the precipitation process does not progress beyond this point in 
naturally aged specimens* 

Lattice parameter measurements confirmed that the above struc- 
tures were present* A constant value of parameter is characteristic 

12 

of the G*P. zones. It was observed that the parameter remained 
constant during natural ageing. The parameter also remained constant 
through two days ageing at 165^ C, indicating that G. p.ta.1 was present 
during this time. At five days ageing the parameter had relaxed to 
the pure aluminum value. Hence 9 became the predominant intermediate 
at this time. 

^Vhen comparing the Patterson origin peaks, Appendix I, Figs. 4 
through 11, it is observed that the as quenched and naturally aged 
specimens have peaks exactly at the origin while those of the arti- 
ficially aged specimens are shifted slightly. Also the origin peak 
for the naturally aged specimen is roughly twice the magnitude of 
all the others. It has been shown that the Patterson peak height 
is not a true measure of the product of the number of electrons, 

> contained in the two atoms separated by the vector that corres- 
1 3 

ponds to the peak. The peak height is directly proportional to 

but inversely proportional to the sum of the cross-sections of 
the electron density peaks for a one-dimensional Patterson peak. 
Therefore a high Patterson peak is not necessarily due to two heavy 
atoms. The height also depends on the product of the electron 

U lbid . 

12 

A. Guinier, Solid State Physics , Vol. 9, p. 341. 

l^Buerger, op . cit . , p. 26. 
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densities of the atoms in the pair. 

The copper atom has a greater electron density than aluminum. 
Copper contains 29 electrons compared to 13 in aluminum and is 11^ 
smaller in atomic size. Hence the presence of copper in an atom 
pair would substantially increase the Patterson peak height for that 
pair. In a specimen fully aged at room temperature a large fraction 
of the copper atoms are assembled in the G.P. zones, which are a few 
atomic layers thick and are on the average less than lOOA apart. 

When looking perpendicular to the zones in the direction 

there would be a certain number of vectors of length **a*’ between 
atom pairs containing one or two copper atoms. There would be a 
smaller number of like vectors in the as quenched specimen, which 
is assumed to be a random substitutional solid solution of copper 
atoms in the aluminum matrix. The number of such vectors would also 
be less in the artificially aged specimens, which initially contained 
and 0 and later only©. This can best be shoim by consid- 
ering the accepted structure of ©, Fig. 2, petge 13. The 6 struc- 
ture is tetragonal with parameters a =* 4.04^ and c = 5.8A, and the 
axes are parallel to the directions. Only one of the 

three possible orientations of G is shown in Fig. 2. When all three 
orientations are considered together, only two produce copper atom 
separations of 4. 04 A in the direction. Where G.P.^ZQ and 

0 occur together the copper atoms are shared between them and the 

o 

total number available to share in atom pairs 4. 04 A apart is less 
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A. Guinier, Solid State Physics , p. 341. 

1 ^ 

•^Silcock, et al ♦ , op cit . , p. 242 
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than in the naturally aged specimen. This is also true when only G 
is present. Hence one would expect zones to produce higher 

Patterson origin peaks than the other intermediate structures. 

As seen in Appendix I the left origin Patterson peaks occur at 
w values slightly less than ’*a^’ for all artificially aged specimens. 

This shift is real and is not due to an error in calculation or 
plotting of the Patterson function. The peak is at .98a for one days ' 
ageing and at .99a for longer ageing times. This is best explained 
by considering the matrix distortion due to each of the intermediate 
precipitates. When a copper atom is substituted for an altiminura 
atom the matrix lattice contracts locally due to the smaller atomic 
size of copper. All three intermediate precipitates are coherent 
with the aluminum matrix. This coherency between matrix and preci- 
pitate lattices produces distortions in the matrix removed from the 
precipitate interface. The extent of the matrix distortion depends 
on the size and shape of the precipitate particle. Thus when G.P.[^l’\ 
zones are formed in their thin disc shape the lattice distortions 
are small and do not extend far from the interface. The average cell 
parameter is then that of the matrix and the Patterson peak remains 
at I'wj = a. 

^Wien G.P.[^Z"\ is present the lattice distortions are more extensive 
than those due to G.P. ea 9 since the G.P.[2^ zones are thicker and have 
larger diameters. The distortions due to © are also extensive because 
of the 0 c parameter of 5.8A. Thus when G.P.ljZi and 0 are present 
together after ageing one day the lattice is distorted by their 
combined effects and the Patterson peak occurs at \^\ = .98a. As 
the G.P.L^^i zones dissolve and 0 grows with longer ageing times the 
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lattice distortion decreases and is due to 0 alone at peak hardness 
of five days. 

In Appendix I, Figs. 12 through 19, are shoim the Patterson 
functions on a large scale. In general, the number of resolved peaks 
increases as the ageing time increases. This is due to the increasing 
number of harmonics used in the synthesis of the Patterson function. 

The number of harmonics used in each case is determined by the point 
at which the coefficients of the broadened curve have decreased to a 
small value and those of the reference curve still have an appreciable 
value. The use of higher harmonics beyond this point results in a 

1 f\ 

greater uncertainty in the coefficients of the corrected curve. The 
intermediate peaks in the above Patterson functions are all quite small 
compared to the origin peaks. The resolution presented by the inter- 

c 

mediate peaks is considerably greater than the 1*5 to 2 A expected 

o 

using x-radiation of 1.542A wave length. Thus the fine detail is 
false and is only a result of the harmonic analysis. No information 
concerning atomic positions can be gained from it. It is observed 
that the fine details of the Patterson functions for five days ageing 
and longer are essentially the same indicating that the..line profiles 
are very nearly the same. 

8. Conclusions and recommendations. 

The scattered intensity distribution in reciprocal space for 

17 

the aluminum-copper alloy is shown in Fig. 3. The scattering 

^^Stokes, op . cit . , p. 389. 

1 7 

A. Guinier, Solid State Physics , p. 339. 
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streaks arotjnd the nodes of the reciprocal lattice are due to the 
destruction of the strictly periodic electron density distribution 
in the crystal* This destruction is caused by the G, P. zones. 




/ 



Fig. 3. Scattering in the reciprocal space of Al-Cu. 

The strealcs pass precisely through the nodes of the matrix and along 
the^^OO^ rows of the reciprocal lattice. The fact that the G.P* 
zones distort the matrix lattice causes the scattering streaks to 
appear only on the high-ongle side of the reciprocal lattice nodes. 

It iras anticipated that the diffuse streaks about the nodes 
trould contribute to the diffraction process in a distinctive iray. It 
should be possible to differentiate from G.P* since the 

length of the diffuse streeks are related to the size of the zone. 
This distinction could not be realized in the present investigation 
as G.P* [2^ zones irere not obtained by themselves on ageing. However 
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it was possible to distinguish GoP.L'3 from© hj observing the 
height and shift of the Patterson origin peak. 

The amount contributed to the diffraction profile by the diffuse 
scattering streaks is uncertain. It is possible that the broadening 
of the profile due to instrumental factors, which were corrected for, 
is greater than that resulting from the diffuse streaks. However, 
since the instrumental broadening remains constant with ageing the 
profile changes due to the diffuse scattering streaks could be detec- 
table. The thermal vibrations of the atoms about their average position 
in the crystal also produce anisotropic diffuse diffraction effects 
about the nodes of the reciprocal lattice. The thermal diffuse 
scattering is related to the elastic properties of the crystal, which 
certainly change during ageing. Hence, the changes in thermal 
scattering should also contribute to the line profile changes during 
ageing. The amount contributed and the best means of detecting it is 
a difficult problem to analyze quantitive ly. There has been no 
attempt in this investigation to do so. 

There are several ways in which the application of this method 

of analysis could be improved. The use of monochromatic x-rays would 

eliminate the profile broadening due to the spectral spread of the 

incident beam. In order to observe the effect of each intermediate 

precipitate on the line profile it is necessary to have them existing 

separately in the alloy. The transition in the aluminum-copper system 

from G.P.t^Q — ► G.P.^Z"\ could best be observed by ageing a 4% copper 

alloy at 110^ C, and from G.P.\^^-*^ 0 by ageing the same alloy at 190^ C. 

The complete intermediate sequence could be achieved with definite 
separations in structures by ageing a 4.5^ copper alloy at 130^0 but 
would require ageing times of about 400 days. 
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APPENDIX I 



PATTERSON PROJECTIONS 

The Patterson functions for the naturally and artificially 
aged specimens are sho^vn on a small scale in Figs. 4 through 11, 
and on a large scale in Figs. 12 through 19. The ordinate is P(w), 
where w is the interatomic vector. The abcissa represents the 
magnitude of w as it ranges from 0 to "a". 
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FIGURE 4« Pbtterson Fimctlen, As quenched specimen 
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FIGURE 5. Patterson Function, aged 65 days at room temperature 
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FIGUBE 6* Patterson Function for specimen aged I day 
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FIGURE 7* Pattorsea Function for apecimen aged. 2 daya 
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FI6UBE 8* Pa.tterson Functl«n for specimen aged 5 day& 
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FI6UBE 9. Faitteraon Functian for specimen aged 9 days 
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FIGOHE 10. Patterson Function for apeciiiLdn ftged 14 days 
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FIGURE II. Pfttterson Function for specimen siged 30 da.y& 
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FIGURE 12:. Patterson Function, asi quenched specimen 
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FIGURE 13. Patterson Function, aged 65 days at room temperature 
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FIGURE 14 » Patterson Function for specimen aiged 1 day 
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FIGURE 15« Patterson Function for specimen aged 2 days 
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FIGURE 16 • F&tterson Function for specimen eiged 5 days 
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FIGDBE 17« Patterson Function for specimen aged 9 days 
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FIGURE 18. Psktterson Fvinctioa for specimen siged 14 days 
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FIGURE 19. Patterson Function for specimen aged SO days 



APPENDIX II 



FOURIER COEFFICIENTS FOR REFERENCE, BROADENED 
AND CORRECTED PROFILES 

The column headings for the tabulated Fourier coefficients 
are as follows: 



N, the harmonic niimber 

ALPHA, coefficient of the real component 
BETA, coefficient of the imaginary component 
INTENSITY for each harmonic is computed by 

k CA.a’' ) 

where L is the number of equal intervals 
between the data points used to describe 
the curve. 

The corrected coefficients are normalized. The arrow in the 
’’corrected** column indicates the number of harmonics used in the 
synthesis of the Patterson function. 
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APPENDIX III 



EXPERBffiNTAL DIFFRACTION LINE PROFILES 



The following instrumental settings were used in obtaining all 
of the diffraction line profiles* The abcissa is in degrees 20* 

Cu Koc radiation, 50 KV, 40 ma, Ni filter 

Scan Speed /minute, Chart Speed 30’’/hour 

Time Constant 2 seconds. Scale 64 x .6 

o 

Angular Aperture 1 , Receiving Slit 0.003" 

Xenon Proportional Counter, 1725 volts 
Base 7 volts, Window 3 volts 
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APPENDIX 



IV 



THE FOLLOWING IS A FORTRAN 60 COMPUTER PROGRAM 
THAT WILL PERFORM THE STOKES UNFOLDING AND PATTERSON 
FUNCTION SYNTHESIS. 



PROGRAM UNFOLD 



C 

C 

C THIS PROGRAM WILL COMPUTE 

C 1. THE FOURIER COEFFICIENTS A(N»1) AND B(N.l) OF A CURVE 

C 2.; THE FOURIER COEFF I C I ENTS A ( N ♦ 2 ) AND B(N,2) OF A SECOND CURVE 

C ' 3.' THE FOURIER COEFFICIENTS A(N) AND 8(N) OF THE UNFOLDED CURVE 2/1 

C A. THE INVERSE FOURIER TRANSFORM (SYNTHESIS) OF THE UNFOLDED CURVE 

C 5. THE CENTER OF GRAVITY OF THE UNFOLDED CURVE 
C ' 6. THE LOCATION OF THE MAXIMUM IN ThE UNFOLDED CURVE 

C 7. A SERIES OF CURVES CAN BE SUCCESIVELY UNFOLDED 

C THE RESULTS CAN BE PRESENTED IN BOTH TABULAR AND GRAPHICAL FORM 

C GIVEN THE COEFFICIENTS I F ( H . K ♦ L ) I **2 ♦ THE PROGRAM COMPUTES THE 
C PATTERSON FUNCTION, Y(U,V,W) (I.E. THE SUMMATION OVER THE COSINE 
C TERMS OF A FOURIER TRANSFORM OF I ( F , ( H , K , L ) I **2 ). 

C THE RESULTS ARE PRESENTED AS ONE OR TWO DIMENSIONAL SECTIONS 
C ONE DIMENSIONAL SECTIONS CAN BE PRESENTED"'! N BOTH GRAPHICAL AND 
C TABULAR FORM 

C TWO DIMENSIONAL SECTIONS ARE PRESENTED IN TABULAR FORM 
C THE EXPERIMENTAL VALUES, I F ( H , K , L ) I **2 ARE CORRECTED FOR LORENTZ, 

C POLARIZATION, MULTIPLICITY AND TEMPERATURE FACTORS 
C 
C 
C 
C 
C 
C 
C 
C 

C • 

C 
C 
C 
C 
C 
C 



FIRST DATA CARD CARRIES GENERAL TITLE - COLUMNS 2-80 
SECOND DATA CARD, 

col' 1-3; NA =' 1 FOR THE FIRST CURVE OF A SET 
- . =2 FOR THE SECOND CURVE OF A SET 
COL 4-5 NB = 0 IF FOLDING IS TO BE CALCULATED 

= i TO COMPUTE FOURIER COEFFICENTS ONLY 
A VALUE FOR NB MUST BE ENTERED OR STOP 15 
WILL RESULT. 

DIMENSIONS OF AA, BB, AND FIX DEPEND ON 
VALUE OF NA SPECIFIED 

COL 6-7 NC = +1 SYNTHESIS OF UNFOLDED DATA 
= -1 SYNTHESIS OF NEW DATA 



5S 



C COL 8-9 NO = 0 NO GRAPH OF PATTERSON FN. 

C =1 GRAPH PATTERSON FN. SEE BELOW. 

C COL 10-74 TITLE OF UNBROADENED CURVE 

C 

C THIRD DATA CARO 

C 
C 
C 
C 

c 

c ^ 
c 
c 
c 
c 

C I 

c 
c 
c 
c 
c 
c 
C 
c 
c 

C FOURTH DATA CARD - UNBROADENED CURVE.N EQUALLY SPACED POINTS FORMAT 

C (16F5.1) , ON THE NECESSARY NUMBER OF CARDS. 

C 

C REPEAT SEQUENCE WITH A DUPLICTE GENERAL TITLE CARD ♦ A SIXTH 

C CARD SETTING NA = 2 AND CARRYING A TITLE FOR THE BROADENED 

C CURVE . A SEVENTH CARO SETTING THE CONSTANTS FOR ANALYSIS OF THE 

C BROADENED CURVE »AND AN EIGHTH CARD WITH THE N POINTS OF THE 

C BROADENED CURVE. 

C 

C FOR A GRAPH OF THE PATTERSON FN, FOLLOW THE TWO SETS OF DATA CARDS 

C (FOR THE REFERENCE AND BROADENED CURVES) WITH TWO MORE CARDS. 

C THE FIRST ONE IS BLANK AND THE SECOND ONE SETS NA=NB=NC=0, ND=1. 

C 

C . FURTHER CONVOLUTIONS CAN BE CARRIED OUT WITH 

C ADDITIONAL SETS OF FOUR DATA CARDS , INCREASING NA BY ONE FOR 

C EACH FOLD UP TO 10. 

C 

C ADDITIONAL CURVES CAN BE UNFOLDED BY ADDING ADDITIONAL SETS OF 

C DATA CARDS SETTING NA *1 AND 2 IN THE SECOND AND SIXTH CARDS OF 

C EACH SET 

C 

C . FOR A FOURIER ANALYSIS SET NB = 1 IN COL 5 IN THE 
C SECOND CARD OF A SET OF DATA CARDS 

C 



- CONSTANTS 
COL 



FOR 



1 




5 


N = 


6 




10 


NCG 


11 


- 


15 


MODE 


16 


- 


20 


J - 1 


21 


- 


25 


K - 1 


26 


— 


30 


L- F 


31 


— 


35 


K2- . 


36 


- 


40 


' LC = 








\ 


41 


— 


50 


THO 


51 




60 


THF 


61 


— 


73 


XL = 



OF UNBROADENED CURVE 
NO OF intervals 

= -1 FOURIER SYNTHESIS - GRAPH 
= 0 NO SYNTHESIS 

= +1 FOURIER SYNTHESIS 
= 1 IMPUT DATA PRINTED 
= 0 IMPUT DATA NOT PRINTED 
FIRST OUTPUT HARMONIC 
(NORMALLY = 0) 

HARMONIC INCREMENT 
(NORMALLY = 1) 

INAL HARMONIC 

SUB HARMONIC - IF 0, PROGRAM 
SETS K2 = 1 . 

0 READS DATA CARD 
MEMORY ■ 

1 USES DATA available IN 
=INITIAL TWO theta' 

= FINAL TWO THETA 






\ 






I 



I 




n n on 



I 



THE PROGRAM IS TERMINATED BY TWO BLANK DATA CARDS 



C 
C 
C 
C 

C STATEMENT NUMBERS 



c 


1 


2 


3 


4 


5 


6 


7 


c 


18 


19 


20 


21 


22 


23 


24 


c 


35 


36 


37 


38 


39 


40 


41 


c 




51 


52 


53 








c 


60 


61 












c 


-- 


\ , 




103 


104 


105 


406 


c 


110 


111 


112 


113 


114 


115 


116 


c 


120 


121 


122 


123 


124 


125 


126 


c 


130 


131 


132 


133 


134 


135 


136 


c 


140 


141 


142 


143 


144 


145 


146 


c 


150 


151 


152 


153 


154 


155 


156 


c 


160 


161 


162 


163 


164 


165 


166 


c 




501 


502 


503 


504 


505 


506 


c 


900 


901 


902 




904 




906 


c 


47 . 


48 49 


' 70 


170 


71 


72 73 



C 

C 

C 



8 


9 


10 


11 


12 


13 


14 


25 


26 


27 


28 


29 


30 


31 


42 


43 


44 


45 


46 







107 


108 


109 


117 


118 


119 


12 7 


128 


129 


137 


138 


139 


147 


148 


149 ■ 


157 


158 


159 


167 






507 


508 


509 510 511 512 513 


907 


908 


'^09 






\ 



15 16 

32 33 

/ 

/ 



514 515 



17 

34 



1 FORMAT (IHl////) 

2 FORMAT (80H0 

1 ) 

24 FORMAT (I5» 2F10.4» F10.5 » 2F10.4 ,F10.5 ♦2F10.4* F10.5) 

39 FORMAT! 115, 2F15.4, F15.4//) 

40 FORMAT! 1H013X1HN10X5„HALPHA11X4HBETA6X 9HINTENSITY) 

900 FORMAT ! 13. 312, 0A8) 

901 FORMAT!8I5, 3F10.3) • 

902 FORMAT !16F5.1) 

904 FORMAT! 3H0N=I4, 7H J=I5. 7H K=I5. 7H L=I5. 

, 1 8H K2=I5.7X6HNCG = I3///12H DATA POINTS//) 

906 FORMAT ! IHO 16X 9HREFERENCE 20X IIHBROADENED 21X 9HCORRECTED //) 

907 FORMAT! IHO 4X IHN 2X 5HALPHA 7X 4HBETA 3X 9HINTENSITY 4X 5HALPHA 
1 ,6X 4HBETA 2X 9HINTENSITY 4X 5HALPHA 6X 4HBETA 2X 9H I NTENS I TY// ) 

908 FORMAT! IX) 

909 FORMAT! IX, 16F5. 1) 



DIMENSION Y!200), AA!200,4), BB!200,4), FIX!200,4), NAME!8), 

1 FR!200,4), FI!200,4), FIT!200,4), 

2 Z!l, 1,900) 




ST 



DIMENSION WW(900). ITIYlEIIZ) 

DO 71 1=1, 12 
■ 71 ITI TLE( I ) = 8H 

ITITLE(l) = 8H LARSEN 
ITITLE(2) = 8H 0278 
ITITLE(7) = 8HPATTERSO 
ITITLE(8) = 8HN FUNCT 
ITITLE(9) = 8HION FOR 
ITITLE(IO) = 8H AL-CU 
LABEL = 4H 65D 

TPI = 6.28318530718 
R=l. 0/(57. 2957795 13 1*2.0) 

3 READ 2 

READ 900, NA, NB, NC , ND, NAME 
IF(NC) 122, 15, 47 
i 47 IF (NA) 15, 15, 9 
I 15 IF (ND) 48,48,49 

49 CALL DRAW( 101, Z,WW, 0,0, LABEL, I TITLE, 0,1 0.0, 0,0, 0,0,9,10.,!, LAST) 
GO TO 3 
. 48 STOP 48 

9 READ 901, N, NCG, MODE, J, K, L, K2, LC , THO, THF, XL 

IF(K2) 30,31,30 i 

; 31 K2 = 1 
30 FN = N 

TNR = 2./FN r 

THF = THF*R 
THO = THO * R 

IF LC EQUALS ZERO, READ NEW DATA, OTHERWISE OPERATE ON OLD 
IF (LC) 32,33,32 

33 READ 902, (Y( I ) , 1= 1,N ) 

32 PRINT 1 
PRINT 2 

PRINT 900, NA, NB, NC, ND, NAME 
IF (MODE) 12,13,12 

12 PRINT 904, N, J, K, L, K2 , NCG 
PRINT 909, (Y( I ) , I = 1 ,N) 

JF (NA - 1) 515i 513, 514 
.515 STOP 515 

513 AREA = 0.0 , 

DO 511 I = 1,N 

;AREA1 = (Y( I )*(SINF(THF) - SINF(THO)))/ (XL * FN ) 

511 AREA = AREA + AREAl 

514 DO 512 I = 1,N 

512 Y( I ) = Y( I ) / AREA 

13 XK2 = K2 
XJ = J 

ARG = TPI/FN/XK2 
ARGl = ARG*XJ 



501 

502 

503 

504 

505 

506 

507 
5,08 

509 

510 
14 
18 



16 

17 



19 



cs = 


COSF(ARGl) 


TCS = 


^ 2.*CS 


SS = 


SINF( ARGl ) 


XK = 


K 


ARG2 


= ARG*XK 


CB = 


COSF( ARG2 ) 


s!b = 

M1 = J 


SINF(ARG2) 



GO TO { 501»502,503»504, 505, 506»507»508, 509*510) 

NAN=NA 

GO TO 14 

NCW=NA 

NAN=NA-1 

GO TO 14 

N;CW = NA-1 

NAN=NA-2 

GO TO 14 

NCW=NA-2 

NAN=NA-3 

GO TO 14 ^ . 

NCW=NA-3 
NAN=NA-4 
GO TO 14 
NGW=NA-4 
NAN=NA-5 
GO TO 14 

NCW=NA-5 ■ 

NAN=NA-6 
GO TO 14 
NCW=NA-6 
NAN=NA-7 
GO TO 14 
NCW=NA-7 
NAN=NA-8 
GO TO 14 
NCW=NA-8 
NAN=NA-9 

BEGINNING OF LOO^ ON Ml 
T3 = 0. 

T2 = 0. 

I = N 
T1 = TCS- 
I = I-l 

IF (I) 16,19,17 
STOP 16 
T3 = T2 
T2 = T1 
GO TO 18 
M=M1+1 



NA 



* T2 - T3 + Y ( I ) 



S<\ 



I 



+ BB(M,NA) * BB(M»NA) )/TNR 



AA(M»NA) = (Tl'- CS^^T2)*TNR 
BB(M»NA) = (SS* T2)*TNR 
FIX(M,NA)=( AA(M»NA)*AA (M»NA) 

TEST FOR END OF (Ml SERIES 
rF(Ml-L) 21, 26,20 

20 STOP 20 

21 M1=M1+K 

T.CS = CB*CS - SB*SS ' 

SS = CB*SS + SB*CS 
CS = TCS 
TCS = TCS + CS 
GO TO 14 

26 NO=NA-l \ 

11 IF(NO) 23,42, 27 

23 STOP 23 

42 IF(NB) 27,3,27 

27 PRINT 1 
PRINT 2 
L2=L+1 
J2=J+1 

IF(NB) 37,38,37 

37 PRINT 40 
PRINT 908 

DO 43 M = J2, L2 
M2=M-1 

43 PRINT 39,M2,AA( M,NA) ,BB (M,NA) ,FIX(M,NA) 

GO TO 3 

38 PRINT 906 
PRINT 907 
IF(NA-2) 22,22,28 

28 AA(M,NAN)=AA(M,NO) 

BB(M,NAN)=BB(M,NO) 

FIX(M,NAN) = ( AA(M,NAN )*AA(M,NAN)+ BB ( M , NAN ) *BB ( M , NAN ) ) / TNR 

22 DO 25 M=J2,L2 
M2=M-1 

J2 = J + 1 
FR(M,NO) = 

1 *TNR 
FI(M,NO) = 

1 *TNR 

IF(M -J2)46, 45, 

45 FA = FR(J2,NO) 

46 FR(M,NO) = FR(M,NO)/FA 
FI (M,NO) = FI (M,NO)/FA 

FIT(M,NO) =(FR(M,NO)*FR(M,NO) + FI(M,NO)^f F I ( M , NO ) ) /TNR 
PRINT 24 ,M2,AA(M,NAN) ,BB(M,NAN)‘,FIX(M,NAN) , AA(M,NCW) 

1) , FIX(M,NCW) ,FR (M,NO) , FI(M,NO) ,FIT(M,NO) 

IF(XMODF(M ,45) )29,35,29 
35 PRINT 1 



(AA(M,NCW)* AA(M,NAN)+ BB ( M , NCW ) *BB ( M , NAN ) ) /F I X ( M ,NAN ) 

f 

(BB(M,NCW)* AA(M,NAN)- AA(M,NCW)* BB ( M , NAN ) ) /F I X ( M , NAN ) 



46 



BB ( M,NCW 



GO 



PRINT 2 
PRINT 906 I 
PRINT 907 
GO TO 25 

29 IF(XMODF(M ,5))25.36.25 
36 PRINT 908 

GO TO 25 

25 CONTINUE V 

IF(NCG) 122, 44, 122 / 

C 

C 

C . PROGRAM FOUSYN 

C 

C 

C 1ST DATA CARD 

C , iCOL 1-80 TITLE 

C : ' 

C 2ND DATA CARD 

C COL 1-4 DEBYE TEMP. B 

C , B=0.0 IF TEMP. FACTOR IS NOT COMPUTED 

C ' 5-11 LAMBDA 

C 11-14 NUMBER OF PEAKS USED 

C COL 15 - 20 NUMBER OF INTERVALS/CELL IN X DIRECTIONS 

C 21-25 NUMBER OF INTERVALS/CELL IN Y DIRECTION 

C , 26 - 30 NUMBER OF INTERVALS/CELL IN Z DIRECTION 

C 

C ADJUST DIMENSION STATEMENT . FOR FRACTIONS UVW OF XYZ USED 

C I 

C , 3RD DATA CARD 
C COL 1-3 MAX VALUE OF U 

C 4-6 MAX VALUE OF V 

C 7-9 MAX VALUE OD W 



c 


10 - 12 


MAX 


value 


OF 


H 






c 


13 - 15 


MAX 


VALUE 


OF 


K 






c ; 


; 16-18 


MAX 


VALUE 


OF 


L 






c ‘ 


(• 19 - 21 


MAX 


VALUE 


OF 


U 


FOR 


PRINTOUT 


c 


21 - 23 


MAX 


VALUE 


OF 


V 


FOR 


PRINTOUT- 


c 


24 - 26 


MAX 


value 


OF 


w 


FOR 


PRINTOUT ACCORDING TO SECTIONS 


c 


WANTED 















C , ■ " 

C ADJUST DIMENSION STATEMENT FOR MAX VALUES OF HKL USED 

C 

C 4TH DATA CARD 

C COL 1 - 3 H 

C 4 - 6 K 

C - 7 - 9 L 

C , \ 10 - 14 RELATIVE INTENSITY (WITH DECIMAL POINT) 

C 15-18 TWO THETA 

C 19-20 MULTIPLICITY 






\ 



c 

c 

122 CONTINUE 
C 
C 

51 FORMAT (3(313* F5.0, F4.0» 12* F5.0)) 

52 FORMAT(F4.2,F7.5, 14, 315 ) 

53 FORMAT (613, 3F4.2, 2F10.4) , 

4 FORMAT (1H0,2X, IHN, 6X, 22H( U , V ’ , W )9X, '9H Y(U,V,W) 

1 // ) 

5 FORMAT (15, 6X, F6.4 , IX, F6.4, IX, F6.4, 6X, F12.4) 

6 FORMAT! 50H ) 

. 7 FORMAT (IX) 

8 FORMAT (IHl//// ) 

60 FORMAT! 21HOCENTER OF GRAVITY = F6.2) 

61 FORMAT (55H0MAX HARMONIC USED FOR SYNTHESIS OF CORRECTED CURVE = 

1' 13 //) 

C 

167 DIMENSION F(l, 1,500), TTH ( 1 , 1 , 50'0 ) , MU(1 , 1,500), C3(100.), 

1 S3(100), FKl, 1,500), Cl(l), C2(l), SKI), S2 ( 1 ) 

C ' 

IF(NC) 129, 166, 128 
166 STOP 166 
129 READ 6 

READ 52 ,B,XL,NP,NU,NV,NW 

READ 53, Jl» J2, J3» UMAX, I2MAX, U3MAX, U1 , VI, Wl, TTHO, TTHF 
GO TO 136 
128 NNl =0 

• ' NN2 =0 ; 

NU = 0 

NV = 0 

NW = N , ' . 

J1 = 0 

J2 = 0 

J3 = 100 

Jll = J1 +1 

J22 = J2 +1 , f 

J33 = J3 +1 ' 

M2 = M - 1 

UMAX = 0 
•I2MAX = 0 
I3MAX = M2 
U1 = 0.0 
V'l = 0.0 
Wl = FN 
TTHO = THO 
TTHF = THF 
136 IlM =' I IMAX + 1 
I2M=I2MAX+1 



ez 



\ 

\ 



n n 



I3M = I3MAX +1 
TTHD = TTHO - TTHF 
DO 115 13=1, I3M 
IS3 = ( 13 - 1 )*100 
DO 115 12 = 1, I2M 

IS2 = ( 12-1 )*10 + IS3 
DO 115 11 = 1, IlM 

IS=IS2+I1 / 

F(IS)=0.0 / 

FI ( IS) = 0.0 
' TTH(IS)=0.0 

MU ( I S ) = 0 . 0 
115 Z( IS)=0.0 

IF (NO 127, 15, 126 

126 DO 137 M = 1, L2 
IF ,(M-1) 159,159,160 

159 I4m' =0 

160 IF(I4M) 161,161,1.58 

161 F( 1 , 1 , M) = FR (M, NO) 

FI Cl , 1, M) = FI (M, NO) 

IF(FR(M, NO) ) 170, 170, 137 

170' I4M = M 
GO TO 137 

158 F(l, 1, M) = FR(M, NO) 

FKl, 1, M) = FKM, NO) 

137 CONTINUE 

FOR TRY = -1.0, Y(U,V,W) IS SYNTHESIZED USING THE FIRST M 
HARMONICS THAT KEEP Y(U,V,W) POSITIVE." 

162 TRY = -1.0 
I3M = I4M -1 
I5M = I4M -2 
GO TO 125 

C FOR TRY = 1.0, Y(U,V,W) IS SYNTHESIZED USING ALL L HARMONICS. 

163 TRY = 1.0 
I5M = L2 - 1 
I3M = L2 

GO TO 125 ■ , • ■ 

127 IF (XMODF(NP, 3)) 111, 112, 111 

112 NP = NP/3 
GO TO 113 

111 NP = (NP/3) +1 

113 DO 114 NR = 1, NP 

READ 51, NNl, NN2, NN3, F(IS), TTH(IS), MU(IS)» FKIS) 

IS=NN3*100+NN2*10+NN1+1 

TTH( IS) =TTH( IS)*R 

AAA=SINF( TTH( IS ) ) 

AAA2=AAA**2 

TE=1.0/EXPF( B#AAA2/XL**2 ) 

, ZZ=COSF(2.0*TTH( IS) ) 
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POL=(1.0+ZZ**2)/AAA2*COSF(TTH( IS) ) 
AMU( IS) = MU( IS) 

FKIS) = F1(IS)/(TE POL *AMU(IS)) 
llA F(IS) = F(IS)/(TE *POL *AMU(IS)) 
125 PRINT 8 

IF (NO 131» 15» 130 
131 PRINT 6 
60 TO 132 
130 PRINT 2 

PRINT 61. I5M 
■132 PRINT A 

TPI=6. 28318530718 
MARK=0 

IF(NC) 157, 138, 139 

138 STOP 138 

139 JU = 1 
UV = 1 

GO TO lAO 

157 DO lOA JU = 1, Jll ' 

DO lOA JV = 1, J22 ^ 

lAO DO lOA JW = 1, J33 
IF(NC)153, 15A, 155 
15A STOP 15A 

155 U = 0.0 

V = 0,0 
60 TO 156 

153 AJU = JU - 1 
ANU = NU 
U = AJU/ANU 
AJV = JV - 1 
ANV = NV 

V = AJV/ANV 

156 AJW = JW - 1 
ANW = J3 

,W = AJW/ANW 
IF(NC) lAl, 1A2, 1A3 
1A2 STOP 1A2 
1A3 II = 1 
12 = 1 
GO TO lAA 

lAl DO 103 11=1, IlM 
DO 103 12=1, I2M 
lAA DO 103 13=1, I3M 



NNl 


= 


I 1 - 


1 


NN2 


= 


12“ 


1 


NN3 


= 


13 - 


1 


FNl 


= 


NNl 




FN2 


= 


NN2 




FN3 


= 


NN3 
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IF (NO 145, 146, 147 

146 S/OP 146 

147 SI ( ID = 1.0 
S2( 12) = 1.0 
GO TO 148 

145 SKID = SINF(U•^^FN1^^TPI ) 

S2 ( 12) = SINF(V*FN2«-TPI ) 

148 S3(I3) = SINF(W*FN3*TPI ) z 

TS = SKID * S2(I2) * S3(I3) * FKI1,I2,I3) / 

IF (NO 149 , 150, 151 

-150 STOP 150 

151 CK ID = 1.0 

C2 ( 12) = 1.0 

GO TO 152 

149 CKID = C0SF(U*FN1<^TPI ) 

C2(I2) = C0SF(V*FN2*TPI ) 

152 C3 ( I3)=COSF(W*FN3*TPI ) 

t = CKID*C2(I2)*C3(I3)*F( 11,12,13) 

IF ( 13 -1 ) 117, 116,^ 117 ' 

116 TSM=0.0 

117 TSM = TSM + T + TS 

103 Z( JU,JV,JW) =TSM 

107 |MARK = MARK +1 
j'jW = JW - 1 

PRINT 5, JJW, U,V,W, Z(JU, JV,JW) 

WW(JW) = JJW 

IF (XMODF (MARK, 35)) 109,108,109 

108 IF (NO 134, 15, 133 

134 PRINT 8 

PRINT 6 
GO TO 135 

133 PRINT 8 

PRINT 2 

135 PRINT 4 "" ... 

GO TO 165 p 

109 IF (XMODF (MARK , 5)) 165,110,165 

110 PRINT 1 . 

1-65 CONTINUE 

104 CONTINUE 

IF (TRY) 72,44,44 
44 GO TO 3 

72 'IF (NO) 163,163,73 

73 CALL DRAW( 1 0 1 , Z , WW , 0 , 0 , LABEL , ITITLE, 0,10., 0,0, 0,0, 6, 10,1, LAST) 
GO TO 163 

END 

END 
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